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The dynamics of a quantum system undergoing frequent "measurements," leading to 
the so-called quantum Zeno effect, is examined on the basis of a neutron-spin experiment 
recently proposed for its demonstration. When the spatial degrees of freedom are duely 
\ taken into account, neutron-reflection effects become very important and may lead to an 

^ ' evolution which is totally different from the ideal case. 
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1 Introduction 



A quantum system, prepared in a state that does not belong to an eigenvalue of the 
total Hamiltonian, starts to evolve quadratically in time [0, §]. This characteristic 
behavior leads to the so-called quantum Zeno phenomenon, namely the possibility 
of slowing down the temporal evolution (eventually hindering transitions to states 
different from the initial one) ||^. 

The original proposals that aimed at verifying this effect involved unstable systems 
and were not amenable to experimental test [Q. However, the remarkable idea to 
use a two- level system motivated an interesting experimental test 0, revitalizing a 
debate on the physical meaning of this phenomenon [0, |[. There seem to be a certain 
consensus, nowadays, that the quantum Zeno effect (QZE) can be given a dynamical 
explanation, involving only an explicit Hamiltonian dynamics. 

It is worth emphasizing that the discussion of the last few years mostly stemmed 
from experimental considerations, related to the practical possibility of performing 
experimental tests. Some examples are the interesting issue of "interaction-free" 
measurements and the neutron-spin tests of the QZE P, [l^]- In practical cases. 



one cannot neglect the presence of losses and imperfections, which obviously conspire 
against an almost-ideal experimental realization, more so when the total number of 
"measurements" increases above certain theoretical limits. 

The aim of the present paper is to investigate an interesting (and often overlooked) 
feature of what we might call the quantum Zeno dynamics. We shall see that a series 
of "measurements" (von Neumann's projections |TT[]) does not necessarily hinder the 
evolution of the quantum system. On the contrary, the system can evolve away from 
its initial state, provided it remains in the subspace defined by the "measurement" 
itself. This interesting feature is readily understandable in terms of rigorous theorems 
[§], but it seems to us that it is worth clarifying it by analyzing interesting physical 
examples. We shall therefore focus our attention on an experiment involving neutron 
spin 1^ and shall see that in fact this enables us to kill two birds with one stone: 
not only the state of the neutron undergoing QZE will change, but it will do so in a 
way that clarifies why reflection effects may play a substantial role in the experiment 
analyzed. 

In the neutron-spin example to be considered, the evolution of the spin state is 
hindered when a series of spectral decompositions (in Wigner's sense [0) is performed 
on the spin state. No "observation" of the spin states, and therefore no projection a 
la von Neumann is required, as far as the different branch waves of the wave function 
cannot interfere after the spectral decomposition. Needless to say, the analysis that 
follows could be performed in terms of a Hamiltonian dynamics, without making 
use of projection operators. However, we shall use in this paper the von Neumann 
technique, which will be found convenient because it sheds light on some remarkable 
aspects of the Zeno phenomenon and helps to pin down the physical implications of 
some mathematical hypotheses with relatively less efforts. 

The paper is organized as follows. We briefly review, in the next section, the 
seminal theorem for the short-time dynamics of quantum systems, proved by Misra 
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and Sudarshan [Q. Its application to the neutron-spin case is discussed in Sec. ^. In 
Sees. I and unlike in previous papers ||, |10|, we shall incorporate the spatial (1- 



dimensional, for simplicity) degrees of freedom of the neutron and represent them by 
an additional quantum number that labels, roughly speaking, the direction of motion 
of the wave packet. A more realistic analysis is presented in Sec. ^. Finally, Sec. ^ is 
devoted to a discussion. Some additional aspects of our analysis are clarified in the 
Appendix. 



2 Misra and Sudarshan's theorem 

Consider a quantum system Q, whose states belong to the Hilbert space Ti and whose 
evolution is described by the unitary operator U{t) = exp{—iHt) , where if is a semi- 
bounded Hamiltonian. Let E he a. projection operator and ETCE = He the subspace 
spanned by its eigenstates. The initial density matrix po of system Q is taken to 
belong to He- If Q is let to follow its "undisturbed" evolution, under the action of 
the Hamiltonian H (i.e., no measurements are performed in order to get informations 
about its quantum state), the final state at time T reads 

p{T) = U{T)poUHT) (2.1) 

and the probability that the system is still in Ti^; at time T is 

P(T) = Tr [u{T)poU^{T)e] . (2.2) 

We call this a "survival probability:" it is in general smaller than 1, since the Hamil- 
tonian H induces transitions out oIHe- We shall say that the quantum systems has 
"survived" if it is found to be in He by means of a suitable measurement process |ll3| . 

Assume that we perform a measurement at time t, in order to check whether Q 
has survived. Such a measurement is formally represented by the projection operator 
E. By definition, 

Po = Ep^E, Tt[poE] = 1. (2.3) 
After the measurement, the state of Q changes into 

Po ^ p{t) = EU{t)poU\t)E/P{t), (2.4) 

where 

P{t) = Tr [u{t)poU^{t)E] (2.5) 

is the probability that the system has survived. [There is, of course, a probability 
1 — P that the system has not survived (i.e., it has made a transition outside He) 
and its state has changed into p'{t) = (1 - E)U{t)poU^{t){l - E)/{1 - P). We 
concentrate henceforth our attention on the measurement outcome ( |2.4| )-( pr^ ).] The 
above is the standard Copenhagen interpretation: The measurement is considered to 
be instantaneous. The "quantum Zeno paradox" |0] is the following. We prepare Q 
in the initial state po at time and perform a series of ii^-observations at times = 
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kT/N {k = 1, ■ ■ ■ , A^). The state of Q after the above-mentioned measurements 

p(^)(r) = Vn{T)poVUt), Vn{T) ^ [EUiT/N)Ef (2.6) 
and the probabihty to find the system in He ("survival probabihty") is given by 

P(^)(T) = Tr K(r)po4(^)l • (2.7) 



Equations (P^)- (|2.71 ) display the "quantum Zeno effect:" repeated observations in 
succession modify the dynamics of the quantum system; under general conditions, if 
is sufficiently large, all transitions outside He are inhibited. 
In order to consider the N ^ oo limit ("continuous observation"), one needs some 
mathematical requirements: define 

V(T) = hm V^iT), (2.8) 

N—KX> 

provided the above limit exists in the strong sense. The final state of Q is then 

p(T) = V(T)poVt(T) (2.9) 
and the probability to find the system in Ti^; is 

V{T) = lim P(^)(T) = Tr \v{T)poV^ (T)] . (2.10) 

One should carefully notice that nothing is said about the final state p(T), which 
depends on the characteristics of the model investigated and on the very measurement 
performed (i.e. on the projection operator E, which enters in the definition of Vn)- 
Misra and Sudarshan assumed, on physical grounds, the strong continuity of V{t): 

lim V(t) = E (2.11) 

and proved that under general conditions the operators V(T) (exist for all real T and) 
form a semigroup labeled by the time parameter T. Moreover, V^(T) = V(— T), so 
that Vt(T)V(T) = E. This implies, by (|]3D, that 

V{T) = Tr [poVt(T)V(T)] = Tr [poE] = 1. (2.12) 

If the particle is "continuously" observed, in order to check whether it has survived 
inside He , it will never make a transition to H — He- This is the "quantum Zeno 
paradox." 

An important remark is now in order: the theorem just summarized does not 
state that the system remains in its initial state, after the series of very frequent 
measurements. Rather, the system is left in the subspace He, instead of evolving 
"naturally" in the total Hilbert space H. This subtle point, implied by Eqs. ( |2.9| )- 
(p.l2|) , is often not duely stressed in the literature. 

Notice also the conceptual gap between Eqs. (p^Tl) and ( p.lO|) : To perform an 
experiment with N finite is only a practical problem, from the physical point of view. 
On the other hand, the iV — oo case is physically unattainable, and is rather to be 
regarded as a mathematical limit (although a very interesting one). In this paper, 
we shall not be concerned with this problem (thoroughly investigated in ||10[) and 
shall consider the N ^ oo limit for simplicity. This will make the analysis more 
transparent. 
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Figure 1: (a) Evolution of the neutron spin under the action of a magnetic field. An 
emitter sends a spin-up neutron through several regions where a magnetic field B is 
present. The detector Dq detects a spin-down neutron: No Zeno effect occurs, (b) 
Quantum Zeno effect: the neutron spin is "monitored" at every step, by selecting and 
detecting the spin-down component. Dq detects a spin- up neutron. 

3 Quantum Zeno effect with neutron spin 

The example we consider is a neutron spin in a magnetic field [Q. (A photon analog 
was first outlined by Peres fl^-) We shall consider two different experiments: Refer to 
Figures 1(a) and 1(b). In the case schematized in Figure 1(a), the neutron interacts 
with several identical regions in which there is a static magnetic field B, oriented 
along the x-direction. We neglect here any losses and assume that the interaction be 
given by the Hamiltonian 

H = fiBai, (3.1) 

/i being the (modulus of the) neutron magnetic moment, and CTj {i = 1, 2, 3) the Pauli 
matrices. We denote the spin states of the neutron along the 2- axis by | t) and | |). 

Let the initial neutron state be po = Ptt = I t)(T I- The interaction with the 
magnetic field provokes a rotation of the spin around the x-direction. After crossing 
the whole setup, the final density matrix reads 

p{T) = e-'^'^poe^''^ = cos' ^^ptT + sm' ^Pii - ^ sino^T(pti - Pit), (3-2) 
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where uj = 2^3 and T is the total time spent in the B field. Notice that the free 
evolution is neglected (and so are refiection effects, wave-packet spreading, etc.). If 
T is chosen so as to satisfy the "matching" condition cosu;T/2 = 0, we obtain 

p{T)=Pii (T={2m + 1)^, meN), (3.3) 

so that the probability that the neutron spin is down at time T is 

P;(T) = 1 (^T = (2m + 1)^, m G . (3.4) 



The above two equations correspond to Eqs. ( p.l| ) and ( |2.2| ). In our example, H is 
such that if the system is initially prepared in the up state, it will evolve to the down 
state after time T. Notice that, within our approximations, the experimental setup 
described in Figure 1(a) is equivalent to the situation where a magnetic field B is 
contained in a single region of space. 

Let us now modify the experiment just described by inserting at every step a 
device able to select and detect one component [say the down (|) one] of the neutron 
spin. This can be accomplished by a magnetic mirror M and a detector D. The former 
acts as a "decomposer," by splitting a neutron wave with indefinite spin (a superposed 
state of up and down spins) into two branch waves each of which is in a definite spin 
state (up or down) along the 2;-axis. The down state is then forwarded to a detector. 



as shown in Figure 1(b). The magnetic mirror yields a spectral decomposition |T2 
with respect to the spin states, and can be compared to the inhomogeneous magnetic 
field in a typical Stern- Gerlach experiment. 

We choose the same initial state for Q as in the previous experiment [Figure 1(a)]. 
The action of M+D is represented by the operator E = [remember that we follow 
the evolution along the horizontal direction, i.e. the direction the spin-up neutron 
travels, in Figure 1(b)], so that if the process is repeated times, hke in Figure 1(b), 
we obtain 

p(^Hr) = r^(T)po4m = (cos^ ffpn = (cos^ ^yp^^^ (3-5) 

where the "matching" condition for T = Nt [see Eq. ( p.3|) ] has been required again. 
The probability that the neutron spin is up at time T, if observations have been 
made at time intervals t {Nt = T), is 

PriT) = (cos^^y. (3.6) 

This discloses the occurrence of a QZE: Indeed, p/^^(T) > pf^~^^(r) for > 2, 
so that the evolution is "slowed down" as increases. Moreover, in the limit of 



infinitely many observations 



piN)^T)--^p{T)=Pn (3.7) 



N 
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and 

V,{T) ^ lim P|^)(r) = 1. (3.8) 

Frequent observations "freeze" the neutron spin in its initial state, by inhibiting 
(A^ > 2) and eventually hindering (A^ oo) transitions to other states. Notice the 
difference from Eqs. ( p.3|) and (|3.4| ): The situation is completely reversed. 



4 The spatial degrees of freedom 

In the analysis of the previous section only the spin degrees of freedom were taken 
into account. No losses were considered, even though their importance was already 
mentioned in p, |T0[ . In spite of such a simplification, the model yields physical insight 
into the Zeno phenomenon, and has the nice advantage of being solvable. 

We shall now consider a more detailed description. The practical realizability of 
this experiment has already been discussed, with particular attention to the N —>■ oo 



limit and various possible losses |]T0| . One source of losses is the occurrence of reflec- 
tions at the boundaries of the interaction region and/or at the spectral decomposition 
step. A careful estimate of such effects would require a dynamical analysis of the mo- 
tion of the neutron wave packet as it crosses the whole interaction region (magnetic- 
field regions followed by field-free regions containing each a magnetic mirror M that 
performs the "measurement"). However, it is not an easy task to include the spatial 
degrees of freedom of the neutron in the analysis; instead, we shall adopt a simplified 
description of the system, which preserves most of the essential features and for which 
an explicit solution can still be obtained. It turns out that the inclusion of the spatial 
degrees of freedom in the evolution of the spin state can result in completely different 
situations from the ideal case, which in turn clarifies the importance of losses in actual 
experiments and, at the same time, sheds new light on the Zeno phenomenon itself. 

Let us now try to incorporate the other degrees of freedom of the neutron state 
in our description. Let our state space be the 4-dimensional Hilbert space Hp ® Hs, 
where Hp = {\R),\L)} and Hs = {\ T)) I i)} are 2-dimensional Hilbert spaces, with 
R{L) representing a particle traveling towards the right (left) direction along the y- 
axis, and j (|) representing spin up (down) along the z-axis. We shall set, in the 
respective Hilbert spaces, 

|->^Q' i->-C)^ i»^Q' <") 

so that, for example, the state |-R |) represents a spin-down particle traveling towards 
the right direction (+?/). Also, for the sake of simplicity, we shall work with vectors, 
rather than density matrices (the extension is straightforward). 

In this extended Hilbert space the first Pauli matrix (Xi acts only on Hs as a spin 
flipper, CTil I) = I I) and ai\ |) = | t), while another first Pauli matrix ri acts only on 
Hp as a direction- reversal operator, ti\R) = \L) and ti\L) = \R). To investigate the 
effects of refiection we assume that the interaction be described by the Hamiltonian 

H = g{l + an){l + /3a,), (4.2) 
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where g, a and /3 are real constants. By varying these parameters and the total inter- 
action time T, the above Hamiltonian can describe various situations in which a neu- 
tron, impinging on a 5-field applied along x-axis, undergoes transmission/reflection 
and/or spin-flip effects. 

It is worth pointing out that the above Hamiltonian incorporates the spatial de- 
grees of freedom in an abstract way: Only the 1-dimensional motion of the neutron, 
represented by L and R, has been taken into account and all other effects (like for 
instance the spread of the wave packet) are neglected. This amounts to consider a 
trivial free Hamiltonian, which can be dropped out from the outset. This may seem 
too drastic an approximation; however, it is not as rough as one may imagine. In 
fact, over the distances involved (a neutron interferometer), the spread of the wave 
packet can always be practically neglected as a first approximation. The introduction 
of the above two degrees of freedom L and R just corresponds to such a situation 
and the simplicity of the model still enables us to obtain explicit solutions for the dy- 
namical evolution. This can be a great advantage. A realization of such a quantum 
Zeno effect experiment is in progress at the pulsed ISIS neutron spallation source. 
Neutrons which are trapped between perfect crystal plates pass on each of their 2000 
trajectories through a flipper device which cause an adjustable spin rotation. Flipped 
neutrons immediately leave the storage system where they can be easily detected (see 
e.g. PI). 



Since the spin flipper ai and the direction-reversal operator Ti commute with each 
other and with the Hamiltonian (|4.2| ), the energy levels of the system governed by 
this Hamiltonian are obviously Et-^^ = g{l + Ta){l + af3) with r, a = ±. Moreover, 
the evolution of the system has the following factorized structure 

If a neutron is initially prepared in state |-R T)) the evolution operator is explicitly 
expressed as 

e-'^^ = t^+ r^Ti + tiai + rmai, (4.4) 

where t^,ti,r^ and are the transmission/reflection coefficients of a neutron, whose 
spin is flipped/not flipped after interacting with a constant magnetic field B, applied 
along the x-direction in a finite region of space (square potential, stationary state 
problem). See Figure 0. These coefficients are connected with the energy levels by 
the following relation : 

By specifying the values of the parameters g,<y,(3 and the total interaction time 
T, one univocally determines t^,ti,r-^ and rj^. Direct physical meaning can therefore 
be attributed to the constants g, a and j3 in ( |4.2| ) by comparison with the transmis- 
sion/refiection coefficients. For example, in order to mimic a realistic experimental 
setup with given values of )f:||,r||, it is enough to obtain the values of g,a and /? 
from ( [4.5|) and plug them in the Hamiltonian ([4.2|) . The model could in principle 
be further improved by making the constant g energy-dependent. We will consider a 
more realistic Hamiltonian in Sec. |6|. 
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Figure 2: Transmission and reflection coefficients for a neutron initiaUy prepared in 
the \R \) state. 
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5 The ideal case of complete transmission 



In the following discussions we always assume that our initial state is \R t), i.e., a 
right-going spin-up neutron, and consider, for definiteness, the case of total trans- 
mission with spin flipped, i.e., = 1, when no measurements are performed. Of 
course, this has to be considered as an idealized situation, since a spin rotation can 
only take place when there is an interaction potential (proportional to the intensity 
of the magnetic field) which necessarily produces reflection effects (with the only ex- 
ception of plane waves). Stated differently, when the spatial degrees of freedom are 
taken into account in the scattering problem off a spin- flipping potential, complete 
transmission is impossible to achieve: There are always reflected waves. Our model 
Hamiltonian ( [4.2[ ) must therefore be regarded as a simple caricature of the physical 
system we are analyzing. Wave packet effects will be discussed in Sec. 

To obtain a total transmission with spin flipped, the evolution operator should 
have the form e~*^^ oc ai, which is equivalent to either 



or 



That is, 



or 



Case i) cosa^^T = sin/5(7T = cosa/^^^T = 0, (5.3) 
Case ii) sin agT = cos PgT = sin aPgT = 0. (5.4) 



(All other cases, such as total reflection with/without spin-flip can be analyzed in a 
similar way.) In both cases, the evolution is readily computed: 

e-'"^\R t) = phase factor x\Rl). (5.5) 

The boundary conditions are such that the neutron is transmitted and its spin flipped 



with unit probability. For the experimental realization, see [T^ . This is the situation 
outlined in Figure |l|(a). 

We shall now focus on some interesting cases, which illustrate some definite aspects 
of the QZE. Let us see, in particular, how the evolution of the quantum state of 
the neutron is modified by choosing different projectors (corresponding to different 
"measurements" ) . 



5.1 Direction-insensitive spin measurement 

We perform now a series of measurements, in order to check whether the neutron 
spin is up. Let us call this type of measurement a "direction-insensitive spin measure- 
ment," for reasons that will become clear later. Refer to Figure 0(a). The projection 
operator corresponding to this measurement is 

E^ = 1-\R 1){R II - |L i I = ^(1 + as), (5.6) 
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Figure 3: (a) Direction-insensitive spin measurement, (b) Direction-sensitive spin 
measurement. 
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that is, the spin-down components are projected out regardless of the direction of 
propagation of the neutron. In this case, after frequent measurements Ei performed 
at time intervals T/N, the evolution operator in Eq. ( p.6|) reads 

Vj^iT) = (E.e-^^^/'^Ey = E,{t, + r,r{)\ (5.7) 



where t| ~ 1 — igT/N and ~ —iagT/N for large [see Eq. ( [4.4|) ]. Taking the 



limit, one obtains the following expression for the QZE evolution operator defined in 
Eq. 

V(T) = lim Vn{T) = e-'^^Eie-"'^^^\ (5.8) 

Interesting physical situations can now be investigated. Choose, for instance, 
gT = n,a = —1/2,(3 = —1, which belongs to Case i) in Eq. ( |5.3| ) [so that, without 
measurements, the neutron is totally transmitted with its spin flipped, as shown in 
Eq. ( |5.5|) ]. When the direction-insensitive measurements are continuously performed, 
the QZE evolution is V(T) = —iEiTi and the final state is 

V{T)\R]) = -t\L]), (5.9) 

i.e., the neutron spin is not flipped, but the neutron itself is totally reflected. This 
clearly shows that reflection "losses" can be very important; as a matter of fact, 
reflection effects dominate, in this example. Notice that this is always an example of 
QZE: The projection operator Ei in ( p.6| ) prevents the spin from flipping. The point 
here is, however, that Ei is not "tailored" so as to prevent the wave function from 
being reflected! 



5.2 Another particular case: seminal model 

Let us now focus on a model corresponding to Case ii) in Eq. (|5.4| ). The choice 
of parameters, e.g. gT = 7T/2,a = 2n,j3 = —1, obviously fulfills these conditions 
for arbitrary integer n. Total transmission with spin flipped occurs again when no 
measurement is performed. 

When direction-insensitive spin measurements, described by projections Ei, are 
performed at time intervals T/N and in the N ^ oo limit, the QZE evolution operator 
in Eq. ( ^.8|) becomes simply Vi(T) = —i{—l)"Ei and the final state is 

Vi(r)|i?T) = -^(-i)"|i?T), (5.10) 

so that the "usual" QZE is obtained. When n = this is our seminal model 
reviewed in Sec. ^. Obviously, the case n = is not rich enough to yield information 
about reflection effects. In the following subsection the case of nonzero n will be 
discussed. 
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5.3 Direction-sensitive spin measurements 

We now consider a different type of spin measurement. Let the measurement be 
characterized by the following projection operator 

E, = l-\Ri){Ri\, (5.11) 

which projects out those neutrons that are transmitted with their spin flipped. Notice 
that spin-down neutrons that are reflected are not projected out by E2: for this reason 
we call this a "direction-sensitive" spin measurement. Refer to Figure 0(b). Even 
though the action of this projection is not easy to implement experimentally, this 
example clearly illustrates some interesting issues related to the Misra-Sudarshan 
theorem. We shall see that the action of the projector E2 will yield a very interesting 
result. For large A^, the evolution is given by 

V2MT) = {E.e-^^^/'^E^f = e-'^^ (l - i^^Z^^^ + 0(l/iV), (5.12) 

where Z = E2{H/g — 1)E2. The QZE evolution is given by the limit 

V2(T) = lim V2.NiT) = e-''^e-''^^E2. (5.13) 

N-*oo 

To compute its effect on the initial state \R t), we note that, when acting on states 
\R I), \L 1) and \L |), which span the "survival" subspace, the Z operator behaves 
as 

/|/?T)\ / a ap\ /|i?T)\ 
Z |L T) = a /? |L T) . (5.14) 
\\Li)) \a(3 ) \ \Li)) 

Let us choose for definiteness (3 = —1, so that 

{Z -l/2f\R^) = e^\R]), (5.15) 
with 9 = v^8aM-T/2. Thus the final state can be readily obtained 



cos gTe+- sin gTO) \R T) + - sin gTe{\L i) - \L ]) 



(5.- 



6) 



Therefore, for a continuous direction-sensitive (namely, E2) measurement, the prob- 
ability of finding the initial state \R t) is not unity. Part of the wave function will 
be reflected, although the neutron would have been totally transmitted without mea- 
surement [see (^)] or with an ".Ei-measurement" [see ( |5.1(J| )]. 

Clearly, the action of the projector E2 yields a completely different result from 
that of El, in (|5.10|) . This is obvious and easy to understand: the state ( ^.16| ) belongs 



to the subspace of the "survived" states, according to the projection E2. Notice also 
that the probability loss due the measurements is zero, in the limit, because the QZE 
evolution ( p.l3| ) is unitary within the subspace of the "survived" states. 
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Figure 4: Spin-up neutron moving along the +y direction with energy E. The mag- 
netic field points to the +x direction and is zero in the region between y'^ and in 
which the measurements will be made. In these field-free regions the wave functions 
are before measurement and | after the measurement. 



6 A more realistic model 

Let us now introduce a more realistic (albeit static) model. Such a model can be 
shown to be derivable from a Hamiltonian very similar to the one studied in the 
previous sections by a suitable identification of parameters (see Appendix A). The 
effect of refiections in the QZE will now be tackled by directly solving a stationary 
Schrodinger equation, which will be set up as follows. 

Let a neutron with energy E = k"^ /2m and spin up {+z direction), moving along 
the +y direction, impinge on regions of constant magnetic field pointing to the x 
direction, among which there are A^ — 1 field-free regions. The thickness of a single 
piece of magnetic field is a and the field-free region has size h. The configuration is 
shown in Figure ^ Thus we have the one-dimensional scattering problem of a neutron 
off a piecewise constant magnetic field with total thickness D = Na. The stationary 
Schrodinger equation is described by the Hamiltonian 

Hz = ^ + fiBainiy), (6.1) 

where n is the modulus of the neutron magnetic momentum, B the strength of the 
magnetic field and 

n() = f^ for ?/<0, y'r,<y<yn, v'n < V (n = 1,2, . . . , N), 
U for yn-i<y<y',, {n = l,2,...,N), 

with yn = n{a + h) and y'^ = yn~i + b, characterizes the configuration of the magnetic 
field B applied along the x-axis. Refer to Figure |. The incident state of the neutron 
is taken to be \ipin} = ]). Let r|(j) be the refiection amplitude for the spin-up 
(spin-down) component. The wave function for ?/ < is written as 

\^o) = e'''\^) + e-''y[r^\^)+r^\i)]. (6.3) 
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Denoting the transmission amplitudes for spin- up and spin-down as t| and ti, the 
outgoing wave function in the region y > y'j^ reads 

\i^N)^e^'nh\^) + ti\i)]- (6-4) 

Since [ui, Hz] = 0, it is convenient to work with the basis |±) = (| |) ± | J,))/\/2, i.e., 
the eigenstates of ai belonging to eigenvalues ±1. For later use we denote r± — ri±ri 

and t-t = t| ± 

In the field-free region, before the point y = rrin where the nth measurement is 
assumed to take place, y'n < y < the wave function is 

Wn) = E {K,y'^'-''-^ + L'^,ae-''^'-'-^) k), (n = 1, 2, . . . , AT). (6.5) 

(T=± 

On the other hand, in the region after the nth measurement, m„ < y < y„, the wave 
function is 

iV'n) = E (i?n,ae^'=(^-^") + L„,,e-^^(^-^")) (n = 0, 1, . . . , TV). (6.6) 

(T=± 

The relation between the amplitudes of the wave functions iV'n+i) ^iid iV'n) at the 
right- and left-hand sides of the nth potential region is determined by the boundary 
conditions at points y„ and y'^. In fact, we have 



M± ^'^'^ , (6.7) 



where the transfer matrix is given by 

^ _ f cos k±a + icos\irj±smk±a — i sinh r^^ sin /c-i-a \ ffi g") 

^ V i sinh rj± sin k±a cos k±a — i cosh 7]± sin /c-|-a y 

with k± — \/k'^ T 2mnB and k/k± = e^±. This can also be expressed in a concise 
way in terms of the Pauli matrices in this two dimensional space, ri, T2 and T3, as 

M± — e^±(-^+'^i)/^e''^='="'^^e~''='=^''''*"'^^^/^. (6.9) 

We clearly see that the above formula contains all the boundary information at points 
yn and y'^: The first and last factors are the kicks exerted at the boundaries of a single 
piece of constant magnetic field, while the central one represents a free evolution with 
relative energy E ^ ^B. 

In what follows, we shall incorporate the measurement processes performed at 
points rrin as some kind of boundary conditions, connecting the primed and unprimed 
wave functions in the field-free region. 
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6.1 Evolution without any spin measurements 

We first consider the case where there is no measurement at all. This enables us to 
set up the notation and rederive some known results (which will be useful for future 
comparison). In this case the primed and unprimed wave functions must be equal 
Wn) — \'^n) in the field-free region. By virtue of Eq. ( |6.7| ) we obtain 



Rn+l,± \ _ AkbT'i i\j f Rn,± 



Notice the boundary conditions i?o,± = l/"\/2 and Ln^± = together with the 
definitions of transmission amplitude Rn,± = e^^'^^t^j ^2 and reflection amplitude 
Lo,± = r±/\/2. After applying the above equation N times, we obtain the following 
relation 



^^ky. J ^ ([iV]±e''^''-M± - [iV - 1]±J , (6.11) 

where [N]± = (g^ — q^^)/{q± — q'^), with q±,q'^ being the two eigenvalues of the 
transfer matrix e^^^'^^M±, which are determined by 

= cos/c6cos/c±a — cosh?7±sin/c6sin/c-i-a. (6.12) 

When there is only a single piece of magnetic field with length a, i.e., = 1, the 
transmission amplitude of a neutron in the spin state |±) is 

—ika 

ta± = 7 — r (6.13) 

cos K±a — I cosh ?7± sm k±a 

as is well known. From Eq. (|6.11| ), for an arbitrary iV > 1, the transmission amplitude 
of the same neutron passing through a magnetic field with a lattice-like structure as 
depicted in Figure 4 can be written as 

For a neutron in its spin-up state | t), the transmission amplitude with spin unflipped 
is then = (t+ + t_)/2 and that with spin flipped is t| = (t+ — 1_)/2. As a result, for 
a spin-up neutron to go through a constant potential of width = D = Na without 
reflection and with spin flipped, i.e., \ti\ = 1, one should require k±D = n±TT or 

E= ^ + fiB= ^ + -\ 6.15 

with n± two arbitrary integers, their difference n+ — n_ being an odd number. In 
this case of complete transmission = 1, the energy E must be larger than the 
potential fiB. The rest of the analysis above, however, is valid also when the energy 
is less than the potential. 
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Now we consider the case where N tends to infinity and the magnetic field pos- 
sesses a periodic lattice structure. The relation ( |6.1CI|) still holds and in order to 
preserve the translational symmetry along the y axis (that is, to keep the Hamilto- 
nian invariant under a translation of (a+b) along the y-axis), one should have = 1 
owing to the Bloch theorem. Equivalently, the trace of the transfer matrix e'^'''''^^M± 
as given in Eq. ( |6.12| ) should be less than one. This determines the energy band of 
the system: those energies that make the absolute value of this trace greater than 1 
will be forbidden, because for these energies or |5'±|~^ becomes larger than one 
and [N]± tends exponentially to infinity when N approaches infinity. For large A^, 
even if there is no periodical structure, there is always some k that makes this trace 
greater than one (e.g. kb + k±a = In). Therefore the transmission probability will 
tend to zero exponentially when becomes large, even though the energy may be 
very large relative to the potential. This shows that refiection effects in presence of 
a lattice structure are very important; as we shall see, this feature is preserved even 
when projection operators are interspersed in the lattice. 



6.2 Direction-insensitive projections 

We consider now the second situation, when direction-insensitive measurements are 
performed at points m„s. By this kind of measurement, the spin-down components 
are projected out and the spin-up components evolve freely regardless whether the 
neutron is travelling right or left. 

The boundary conditions imposed by this kind of measurement at point m„ for 
the wave function \ipn) and {ipn) in the field-free region are expressed as 

= ^n,i = 0, ( f f 'I ) = e-^^^- ( ^^'1 ) , (6.16) 

where Rn,-\ = {Rn,+ + Rn-)/2 and Rn,i = {Rn,+ — Rn-)/'^ for right-going compo- 
nents and similar expressions for the left-going and primed components. Therefore, 
application of Eq. ( |6.7| ) times yields 

where the 2x2 transfer matrix Mi has the following matrix elements 

(Ml),,- = M,, - AM,2AM2,/M22 (6.18) 

with M = (M+ -I- M_)/2 and AM = (M+ - M_)/2. 

Now we take the limit as required by a "continuous" measurement, i.e., —>■ oo, 
a —>■ keeping Na = D finite and Nb 0. By the definition (|6.8|) of the transfer 
matrix, we have the small-a expansions 

M = l + ikaT^ + 0{a^), AM = Cka{r2 - ir^) + 0{a^) (6.19) 
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with ( = fiB/2E, obtaining 

lim (e^'=''"Wi)^ = e'''^^\ (6.20) 



Recall that = e~'^^^RN,^ is the transmission amplitude, Lo,| = the reflection am- 
plitude and Lat I = and i?o,T = 1 because of the boundary conditions. After taking 
the limit ^ oo in Eq. ( |6.17] ), we see that the transmission (survival) probability 
becomes one, i.e., = 1, for any input energy and magnetic field. This reveals 
another aspect of neutron QZE: When the energy of the neutron is smaller than the 
potential, the transmission probability decays exponentially when the length increases 
and no measurement is performed; By contrast, when continuous direction-insensitive 
measurements are made, one can obtain a total transmission! 

If we choose the energy of the neutron and the potential as in Eq. (|6J[5| ), without 
measurements the neutron will be totally transmitted with its spin flipped. On the 
other hand, if the spin-up state is measured continuously, the neutron will be totally 
transmitted with its spin unfiipped. This is exactly the QZE in the usual sense. Our 
analysis enables us to see that two kinds of QZEs are taking place: One is the QZE 
for the right-going neutron, by which we obtain a total transmission of the right-going 
input state, and another one is for the left-going neutron, which preserves the zero 
amplitude of the left-going input state. This case corresponds to projector Ei in our 



simplified model in Sec. 5.2 



6.3 Direction-sensitive projections 



The third case we consider is the direction-sensitive measurement. By this kind of 
measurement the left-going components (or the reflection parts) evolve freely, no 
matter whether spin is up or down, and the right-going components are projected to 
the spin-up state. The corresponding boundary conditions are 



R 



■n,i 



0, 



-ikb T I 



(6.21) 



If we apply Eq. ( |6.7[ ) times, supplemented with these boundary conditions, the 
following relations among the transmission and reflection amplitudes are obtained 



AkD 




Jkbi: 




(6.22) 



where S3 is a diagonal matrix E3 = 
M2 is given by 

/ Mil 

M2 = 



M12 

M21 M22 
V AM21 AM22 

In the limit of continuous measurements (A^ 00, a 
constant, and Nb 0), the transfer matrix is expanded as 

M2 = 1 - ika/3 + ikaZ2 + O(a^), 



diag{l, —1,-1} and the 3x3 transfer matrix 

(6.23) 

0, while keeping D = Na 
(6.24) 



AM12 
AM22 
M22 
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for small a, with {( = fiB/2E) 



/4/3 -C 
Z2= -2/3 C I , (6.25) 
V C C -2/3, 

and we have 

lim (e^^^^Ws)^ = e-^^^/3gifcZ)Z2_ ^g_26) 

Notice that the matrix Z2 satisfies S3Z2S3 = Zl, from which we obtain, in the 
above limit, the conservation of probability 

i%r+hr+hi' = i- (6-27) 

This shows that there are no losses caused by the continuous direction-sensitive mea- 
surements. On the other hand, the transmission amplitude with spin unflipped is 
explicitly given by 

-i4kD/3 

*T = jpin^^^ (6-28) 

which implies that the transmission probability |t|p is in general not equal to one. 
To have a general impression of its behavior, we plot T| = as a function of kD 
and ( in Figure 

Some comments are in order. There are two critical values for (, namely and 
Cc = 4a/3/9 ~ 0.77. When < C < Cc the matrix Z2 has three real eigenval- 
ues and the transmission probability will oscillate depending on kD. When C = Cc 
the transmission probability will decay according to (kD)^"^. In fact, if one defines 
G = Z2-2/3, it is easy to show that e'"''^^ = 1 - ikDG + {e^'''^ - 1 - 2ikD)G^/A, 
because G satisfies G^{G + 2) = 0. Then one can explicitly confirm that the element 
(e~**'^'^)ii includes a linear kD term, which gives the {kD)~'^ behavior to the trans- 
mission probability. Finally, when C > Cc the matrix Z2 has two imaginary eigenvalues 
and therefore the transmission probability decays exponentially with kD. This can be 
seen clearly in Figure 5(a). An interesting case arises when we consider 1/2 < < (c 
or E < jjiB < 8-\/3-E/9 ~ l.hE. Without measurements, the transmission probabil- 
ity decays exponentially when the length of the magnetic field is increased, because 
the input energy is smaller than the potential. When continuous measurements are 
performed, however, the transmission probability will oscillate as the length of the 
magnetic field increases. 

As we can see in Figure |], although the conditions (|6.15| ) for total transmission in 



absence of measurements have been imposed, the transmission probability T| is not 
one, as it would be for the "ordinary" QZE. Refiections are unavoidable. This case 
corresponds to the projector E2 considered in the simplified model. 

As we have seen, there are peculiar refiection effects in presence of projections, 
when D (total length) is varied. This is clearly an interference effect, which can lead 
to enhancement of reflection "losses," if the "projection" does not suppress the left 
component of the wave (this is what happens for E2). This proves that reflection 
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kD 200 



300 1 




Figure 5: The transmission probability with spin unflipped Tup = is plotted as a 
function of kD and z = C, m. (a) and as a function of Bi = -/mflBD and kD in (b). 
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Figure 6: The transmission probablity with spin unflipped T-^ = as a function 
of n, when the conditions ( |6.15|) for total transmission are satisfied with n_ = n and 

n,=n + 9. 
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effects can become very important in experimental tests of the QZE with neutron 
spin, if, roughly speaking, the total length of the interaction region "resonates" with 
the neutron wavelength. It is interesting that such a resonance effect takes place 
even though the dynamical properties of the system are profoundly modified by the 
projection operators, in the limit of "continuous" measurements, leading to the QZE. 

Finally, we would like to stress again that we are performing an analysis in terms 
of stationary states (i.e., transmission/reflection coefficients for plane waves), while 
at the same time we are analyzing a quantum Zeno phenomenon, which is essentially 
a time-dependent effect. This is meaningful within our approximations, where the 
wave-packet spread is neglected and the measurements are performed with very high 
frequency. A more sophisticated argument in support of this view is given in Ap- 
pendix A. In the present context wave packets effects, if taken into account, would 
result in a sort of average of the effects shown in Figures 5 and 6 (which refer to 
the monochromatic case); however, our general conclusions would be unaltered. It is 
worth stressing that, in neutron optics, effects due to a high sensitivity to fluctua- 
tion phenomena (such as fluctuations of the intensity of the magnetic field) become 
important at high wave number and constitute an experimental challenge p!7[] . 

7 Summary 

We have analyzed some peculiar features of a quantum Zeno-type dynamics by dis- 
cussing the noteworthy example of a neutron spin evolving under the action of a 
magnetic filed in presence of different types of measurements ("projections"). 

The "survival probability" depends on our definition of "surviving," i.e., on the 
choice of the projection operator E. Different Es will yield different final states, and 
Misra and Sudarshan's theorem simply makes sure that the survival probability is 
unity: the final state belongs to the subspace of the survived products. 

In the physical case considered (neutron spin) , our examples clarify that the prac- 
tical details of the experimental procedure by which the neutron spin is "measured" 
are very important. For example, in order to avoid constructive interference effects, 
leading to (unwanted) enhancement of the reflected neutron wave, it is important to 
devise the experimental setup in such a way that reflection effects are suppressed. 
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Appendix A 



In this appendix, we shall endeavor to establish a connection between the models 
analyzed in Sees. ^ and In other words, we will examine whether the parametriza- 
tion of the Hamiltonian of the form ( |4.2| ) is compatible with the more realistic one 
considered in (|0|) and in such a case find which values are to be assigned to the 
parameters and g. To this end, it is enough to consider the scattering (i.e., the 
transmission and reflection) process of a neutron off a single constant magnetic field 
B of width a. We compare the scattering amplitudes calculated on the basis of the 
simple abstract Hamiltonian ( ^.2| ) and of the more realistic one ( |6.1| ). Notice that 
the process is treated as a dynamical one in the former case (T is regarded, roughly 
speaking, as the time necessary for the neutron to go through the potential), while 
in the latter case we treat it as a stationary problem. 

Observe first that the tranfer matrix M± in ( |6.8| ), derived for the stationary scat- 
tering process, yields the following transmission/reflection amplitudes 



2 V(M+)22 (M_) 



22 , 



■(M+)21 ^ (M_)2l' 



2 V(M+)22 ' (M_)22y' ' 2 V(M+)22 

It is easy to show that the relations (|A.1| ) are equivalent to 




(A.l) 



( \ 









\ 


1 -1 


1 -1 






1 1 


-1 -1 






V 1 -1 


-1 1 / 







(A.2) 



where we have introduced 



1 - (M±)21 

(M 



±J22 



(A.3) 



(M±)22 

It is important to realize that these quantities are just phase factors. In fact, since 
(M±)2i = i sinh 'r]± sin k±a and (M-|-)22 = cos k^-a — i cosh ri± sin k±a 

and 



|1 ± (Af±)2i|^ = |(M±)22r = 1 + sinh^ sin^ A;±a, 
their absolute values are unity. Thus we can rewrite them in the form 



M±+<f>±) 



,i(-€±+0±) 



where 



^-t = tan "'^(sinh?7± sin /c-i-a) and 0± = tan ^ (cosh r^^ tan A;-|-a). 



(A.4) 
(A.5) 

(A.6) 
(A.7) 
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Observe now that 
equivalent to 



dynamically derived from the abstract Hamiltonian 



IS 







/I 1 1 
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/e 
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1 -1 1 


-1 
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1 1 -1 
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[l -1 -1 


1 / 


Ve 



-iE+-T 
-iE-+T 
-lE T 



(A.8) 



The apparent similarity between the above relation and ( |A.2| ) , valid in the stationary 
scattering setup, induces us to look for a more definite connection between the two 
cases. 

If we slightly generalize the abstract Hamiltonian (|4.2| ) 



dyn 



g[l + an + (3ai + 7ricri], 



(A.9) 



by introducing the additional parameter 7, we easily find the correspondence existing 
between the parameters involved: The incident wave number k of the neutron and the 
configuration of the static potential (strength B and width a) determine the scattering 
data, which are reproducible by an appropriate choice of parameters a, /3, 7 and gT 
in the dynamical process governed by the Hamiltonian ( |A.9| ). 

For definiteness, consider the case of narrow potential, that is, a — > or fca -C 1. 
Incidentally, notice that this is the case of interest for the analysis of the QZE. The 
above ^± and 0± are then approximated as 



(A.IO) 



where we set C, = fiB/2E = nifiB/k'^, as in Sec. |^. In the limit a — 0, the evolution 
time T is also considered to be of the same order of a and the transmission and 
reflection coefficients are expressed, in terms of the parameters 7 and gT, as 



( 



\ 



1 

iPgT 
iagT 
\-ii9T J 



{All) 



In the stationary scattering problem, the same quantities are calculated to be 

-ika ] 



h 
Vn/ 



ika D/ 



V 



fl-ika + t{(j)+ + (f)^)/2\ 
z(0+-0_)/2 



V 



-z(e+ + ^_)/2 
- e-)/2 



1 

—iQka 


V —iQka J 



(A.12) 



Therefore, the following abstract Hamiltonian 

ifdyn = + ri)(Ti 



(A.13) 



can reproduce the desired scattering data when the system evolves under this Hamil- 
tonian for time T = a/v = ma/k. 
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It is also interesting to see how such a dynamical Hamiltonian ifdyn may reproduce 
the transfer matrix M± ( |6.8|) , which further confirms the equivalence between the two 
formalisms, stationary and dynamical, governed by the Hamiltonians Hz and -ffdyn, 
respectively. For this purpose, consider first a neutron, initially prepared in state 
|i?±), subject to the dynamical evolution engendered by -ffdyn for time T = ma/k. 
By definition, the transfer matrix connects the scattering products in the following 
way 

^\A=mJ^^ ]. (A.14) 



; "vw. 

These scattering amplitudes are given by the corresponding matrix elements of the 
evolution operator e~*^"^, 

g-ifca^/^^ = (i? ± |e-^^'^|i?±), Lo,± = {L± |e-'^'^|i?±), (A.15) 

which reduces, for small T, to 

^ r-^l+ikaT if^BT, Lo,± ~ Tif^BT. (A.16) 

On the other hand, if a neutron is prepared in \L±), we have the relation 

where 

R[ ^ = e'^^'^R ± \e-'"^\L±) ~ T^/i^T, Lq ± = (L ± \e~'"^\L±) ~ 1 t ^f^BT. 

(A.18) 

It is now an easy task to determine the matrix elements of M± from the above relations 
(|AA4D -(|A1^). We obtain 

By defining a "generator" Gd 

Gd = fiB(iT2 + T3)a, - 2Ets, (A.20) 
the transfer matrix M± for finite a (or T) can be rewritten as 

M± = (±|e-^^<^^|±), (A.21) 



which is nothing but the transfer matrix (|6.8| ), obtained for the stationary-state prob- 
lem from the Hamiltonian Hz- 
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